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ABSTRACT 



The major features of a proposed approach to radar bomb 
scoring are discussed. The method uses observed deviations from 
desired release conditions as a basis for predicting bomb mean 
point of impact in the target plane. Circular Error Probable (CEP) is 
estimated using a noncentral chi-square approximation of the bomb 
impact distribution. Contributions of bomb scoring system errors 
(such as errors in radar location) to the apparent delivery 
system CEP are discussed. 



1. Introduction. 



We consider the problem of estimating air-to-ground delivery 
accuracy from radar observed release conditions, rather than impact 
data. Such procedures can relieve requirements of actually dropping 
bombs in training or test situations, with corresponding reduc- 
tion of restrictions on resources such as bomb run location or 
test time required to complete system evaluations. In what follows, 
we present an approach which does not require extensive background 
and training (mathematical, computational, or procedural) on the 
part of bomb scoring personnel. 

Motivation for much of what follows was provided by related 
circular error probable (CEP) and bomb dispersion analyses for the 
AN/TPQ-27 radar bombing system which were conducted at the Naval 
Postgraduate School [1,2]. The proposed procedure suggests use 
of a computer assisted, radar controlled bombing system, with the 
AN/TPQ-27 application serving as a prototype example. The proposed 
method envisions a computer software package for the AN/TPQ-27, 
or similar system, which could provide a readout of release condi- 
tion errors or deviations from predetermined release conditions. 
Through ballistic considerations, these deviations are then trans- 
lated into range and deflection "aim" errors on the ground. Com- 
putation of the estimated CEP is accomplished by means of a model 
based on the non-central chi-square distribution, where the non- 
centrality parameter is a function of the computed range and 



deflection aim errors. 



since in theory, release conditions required to make the 
mean point of impact (MPI) coincident with a given point on the 
ground are not unique, it has been suggested that delivery system 
performance might be evaluated in terms of the envelope of all 
possible combinations of such conditions. This envelope is in a 
space of a dozen or more dimensions, counting components of posi- 
tion, velocity and acceleration; it appears to be difficult to 
determine these envelopes and manipulate them so as to determine 
median miss distance. As we shall see, it is relatively simple 
to relate errors in those components to MPI location in the target 
plane, at least to a reasonable approximation. Since the latter 
approach seems not only more tractible, but more natural as well, 
we have abandoned the "envelope" approach. Keene [7] has discussed 
a process for transforming "small" errors in certain components of 
delivery conditions, to MPI estimates. In what follows we outline 
the principle ideas involved. In addition, we discuss effects of 
errors, such as radar location error, which are outside of the 
delivery system itself. In section 3, we discuss approaches to 
extracting such errors from the CEP estimate, so they would not 
be erroneously charged against the system under test. 

2 . Transforming Delivery Conditions to the Target Plane . 

Except for transient effects such as turbulence, it should 
be possible to measure quite accurately most of the forces acting 
on a bomb at the time of its release. To do so, however, would 
require extensive instrumentation aboard the delivery aircraft 
and a means of transmitting these measurements to the ground for 
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analysis. Thus, the most attractive attributes of radar bomb 
scoring, namely, low cost and flexibility, might be lost in an 
attempt to obtain all of the measures theoretically possible. 

An alternative approach is to base the trajectory estimation 
on data which can be obtained from a radar, or radars, tracking 
the delivery aircraft. This greatly simplifies the computational 
complexity of the MPI estimation problem by limiting the number of 
parameters to be considered in the equations of bomb motion. We 
assume the radar can provide information on the coordinate loca- 
tion of the aircraft at time of release, as well as aircraft 
velocity and acceleration. In addition to tracking information, 
we assume other "inexpensive" inputs are available from sources 
outside the radar-aircraft system. Meteorological data and the 
effects of earth curvature and rotation may be input as standards for 
the location of the bomb MPI. In addition, bomb parameters, such 
as drag curves, may be available for the type bomb being dropped. 

The MPI may be estimated in either of two ways, which we 
call "ballistic computation" and "bombing tables." The first of 
these utilizes ballistic equations to calculate a predicted bomb 
path and its resultant estimated MPI. Ballistic equations of motion 
usually involve consideration of aerodynamic forces, variable winds, 
density and temperature variations, the effects of earth curvature 
and rotation and bomb parameters such as weight, diameter and 
configuration. As mentioned above, these parameter values, neces- 
sary for solution of the ballistic equations, may be available from 
sources outside the bomb scoring system. As is pointed out by 
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McShane, Kelley and Reno [8] , computer assistance is required 
for the numerical solution of such equations . 

The complexity of the ballistic computation model may make 
it too cumbersome for some applications. If so, an alternate 
approach can be used, which takes advantage of the fact that much 
of the work in computing trajectories from models of this type 
has been accomplished and documented in the form of trajectory 
and bombing tables. The use of these tables provides a quick and 
computationally simple means of estimating the MPI when the appro- 
priate corrections for existing local conditions are applied to 
the tabled values. Furthermore, tables are available for nearly 
all ordnance of interest for a bomb scoring system [9] . The 
rationale for suggesting use of these tables in precalculating 
desired bomb release conditions is that this approach is simple 
(inexpensive) and is widely accepted, in that their use is common 
practice in combat situations even where there may be critical 
accuracy requirements. 

With either of the two methods of MPI estimation mentioned 
above (ballistic computation or bombing tables) , evaluations of 
bombing system accuracy can be made for either a pilot controlled 
or computer controlled mode of bombing. In the former case, the 
radar measurements of release conditions provide points of entry 
into the appropriate tables or input values for ballistic computa- 
tions. The tables provide range as a function of release altitude 
above the target and release velocity under assumed atmospheric 
conditions. To these tabled values, corrections due to non-standard 
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conditions may be applied. With the direction of bomb release 
available from tracking data, the estimated point of impact may 
be determined. Comparison of this point with the target location 
yields the desired estimates of range and deflection components 
of the estimated MPI. With the ballistic computation method, an 
MPI is estimated from observed tracking data through a computer 
solution of the ballistic equations. 

It should be noted that for a computer controlled mode of 
bombing, the scoring system could be physically the same as the 
control system. Such a configuration is not as desirable in some 
respects as one using separate physical systems, because with the 
former, control system errors may not be detected. Also, it is 
possible to reduce scoring system errors by choosing an advantageous 
location near the target. Nevertheless, control and scoring with 
a single physical system might be of interest for reasons of economy. 
In the case of computer controlled bombing and scoring with a 
common physical system, the procedure would be somewhat different. 
Local wind conditions and atmospheric data, target and radar data 
and weapon ballistics are preset inputs to the control computer. 

The desired release conditions calculated with these data provide 
and estimate of the release conditions which would place the 
expected point of impact on the target. If it can be assumed that 
deviations from the desired conditions at release are relatively 
small and can be measured to within small error, only the magnitudes 
of these deviations need be considered in estimating range and 
deflection errors. The restriction that these deviations be small 
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is necessary to insure that both the desired and achieved 
trajectories are subjected to very nearly the same conditions and 
forces. The trajectories will then be theoretically nearly 
identical in shape, so the release deviations may be translated 
through simple relationships to range and deflection aim errors 
in the target plane. It would be possible to pre-calculate and 
tabulate range and deflection components of estimated MPI as a 
function of deviations in actual release conditions from desired 
release conditions. 

A detailed discussion of how these deviations are translated 
and combined into MPI estimates is presented and discussed by 
Keene [7] . Basically, the procedure involves elementary geometri- 
cal and positional transformations, and requires the following 
assumptions : 

the desired release conditions are known from the computer 
solution of the ballistics model. 

the desired conditions will place the mean point of impact 
on the target. 

deviations from the desired release conditions are measur- 
able (with small error) and may be output from the bombing 
system. 

deviations are small. 

Since many of the parameters mentioned can be measured 
fairly precisely, the ultimate accuracy of the predicted impact 
point will probably depend on the accuracy with which the radar 
can measure the release conditions: location, velocity, and 

acceleration. The degree of accuracy obtainable is a func- 
tion of the radar being used and will vary from one situation to 
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another. A detailed discussion of radar errors and their 
determination is beyond the scope of this report. In general, 
these errors in measurement are related to the aircraft movement 
relative to the radar as well as the relative positions of the 
radar and aircraft. The use of multiple radars could reduce 
tracking errors, at the expense of operating cost and complica- 
tion of the measurement problem due to the difficulties in cali- 
bration and collimation. In addition, the placement of the tracking 
or "scoring" radar (or radars) relative to the target should be 
selected so as to reduce tracking errors, near the time of release, 
as much as possible. The error characteristics of the radar (or 
radar type) employed may be well known. In a future report we 
shall discuss how knowledge of radar error characteristics can be 
used in an adjustment of predicted CEP, which would avoid "charging" 
radar bomb scoring (tracking system) errors against the delivery 
system. In the present report, it is assumed such errors are small 
relative to delivery system and ballistic errors, so significant 
errors are not caused by ignoring them in the calculated CEP. 

3 . Effects of Errors on Predicted CEP . 

We next discuss how errors, including MPI's offset from 
the actual target location, affect CEP values. We shall then indi- 
cate briefly how estimates of CEP might be adjusted to take such 
errors and offsets into account. In the foregoing section we 
described how the MPI associated with a given single bomb run 
could be estimated from tracking and other data. The amount of 
offset in a given run of the MPI from the target can be thought 
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of as bias in the delivery system, for the given run. We desire 
to relate CEP to bias and ballistic dispersion. Unfortunately, 
there is not a simple relationship (for even a one-dimensional 
version of the problem) , such as the familiar one involving mean 
square error, variance, and bias. 

If there were no differences in range and deflection vari- 
ances, due to ballistic errors, and if these components of ballistic 
error were normally distributed and independent, the ballistic 
distribution would be the so-called "circular normal distribution" 
(if the MPI coincided with the target) or "offset circular normal 
distribution" (if bias was non-zero) . In the first case, squared 
radial error (properly normalized) is chi-square distributed and 
in the second case it is non-central chi-square distributed. In 
mathematical terms, suppose L = (^) is the two-dimensional 
impact location random vector whose first component represents 
deflection miss and whose second component represents range miss. 

It is well known [5] that if L is distributed as described above 

(written "L~N (U/(J^I) ") / then certain quadratic forms in L have 

L *GXj 

non-central chi-square distributions. We write this "—^- 5 — 

^(k X)'" '^here X = y'Gy/2a^ is the noncentrality parameter of 

the non-central chi-square distribution with k degrees of freedom. 

This is the distribution of L'GL/a^ if and only if G is idem- 

potent and of rank k. In particular, with G = I, it follows 

that squared radial miss distance L'lL = normalized to 

L * L 

variance units, — — , is distributed , /- 2 \ • This fact 

^(2,y'y/2a^) 

can be used to compute CEP^ (and hence CEP) , which is the median 
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of this distribution. If y = 0, i.e. MPI is coincident with 
the target, the noncentrality is zero and the distribution of 
L*L/a^ is X( 2 ) ' which coincides with the exponential distribu- 
tion with parameter X = 1/2. Hence, by definition. 



75 - = P[L'L ^ CEP^] 



= P[ 



L'L ^ CEP' 



■] = 1 -e 



- ^EP^/a^ 



which in turn implies the familiar formula 



CEP = /2 Jin 2 0 1.1774 0 . 



( 1 ) 



If y / 0, the median of the non-central chi-square distri- 
bution corresponds to CEP^. Unfortunately, this median cannot be 
expressed in simple closed form as before. We discuss below, in 
section 4, numerical methods of approximating this value. Moreover, 
in practice it is unlikely that the components of ballistic dis- 
persion have equal variances. Suppose L~N 



f 


r^D] 




f^D Oil 


i 




/ 


0 

Q 



0 1 

Then - 2 - + — 2 - is distributed as Y/o ^\/ where \ = rr 

0Q 0 ^ ^(2, A) 2 



'D 



21 






Unfortunately, with 0 ^ / 0 ^, squared radial miss distance, 

Jj H 

, can no longer be normalized to a non-central chi-square 
random variable; rather, various approximations and numerical com- 
putations may be used to approximate the value of CEP in this case, 
as described below in section 4. 

The point of the preceding discussion is that bias (i.e. 

MPI offset from the target) affects the distribution of squared 
radial miss distance through its non-centrality term, at least to 
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a good approximation. We shall refer to the bias (fixed offset y) 
as a fixed effect, and contrast it with random effects, as follows. 
Suppose that for a certain bomb run situation the MPI is y, and 
for a given bomb run in the situation the actual aimpoint is an 
unobservable random value A, where A~N(y,L^). Since A is 
a random variable, aimpoint offset from y is called a random 
effect . Now if a bomb were actually released on this run, the 
location of impact L, would depend on the (unobservable) value 
of A for that run, as well as ballistic error. 

We now give a discussion which arrives at the principle that 
fixed and random effects may affect CEP in different ways, analogous 
to their impact on X. For fixed effects, such as delivery system 
radar location survey error, delivery system radar alignment error 
and drag curve error, the effect may show us as a fixed bias effect, 
that is, an effect upon the mean y. For random effects, that is, 
error components that would theoretically vary randomly from one 
bomb run to another conducted under fixed conditions in a given 
situation, the effects enter through additive variance factors 

Suppose the conditional distribution (ballistic distribution) 
of L for given A is N(A,Zj^), where represents parameters of 

ballistic dispersion in deflection and range; for example, 

in a model described above. The unconditional dis- 

) 

tribution of L, obtained by integration of the conditional density 
with respect to the aimpoint density, is N(y,Zj^+Z^). Thus 

L'GL'\^X(^ X) ^ ^ rank (G) and X = j y'Gy if and only if 

G(Z^+Z,) is idempotent. In particular, with G = (E^+Zj.) 

LA Jj 
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L' (Ej^+E^) , with A = Y y ' ^VJ. While it is the 

squared radial miss distance L'L whose median we seek (we discuss 
this below), rather than that of consideration of 

the latter and its distribution provides insight into how fixed 
and random effects may also affect CEP. Note, for example, that 
the random effect A affects the non-centrality parameter X above 
through the inverse total variance Recall the fixed 

effect, y, affects X through a quadratic function of the compon- 
ents of y. 

The significance of the preceding discussion for the radar 
bombscoring problem is that the procedure used in calculating CEP 
should account for detected foxed effect errors through the esti- 
mated y and for known random effect errors through the accumulated 
variance-covariance matrix E +E +... . Depending on the configu- 

ration of the scoring system relative to the delivery system, some 
of these errors may not contribute to the CEP calculated, even if 
they are detected or are known in advance. For example, suppose 
the scoring system is the same as that used for delivery. Then any 
survey in radar location, radar alignment error, or wind condition 
error will affect delivery and scoring in identically the same ways, 
so nothing would be gained in "correcting" delivery and scoring. 
Rather, in this configuration, only deviation y of measured and 
desired drop conditions are used, together with the ballistic 
variance-covariance, E^ . On the other hand, if the delivery 
system were phvsicallv different from the scoring system, wind 
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condition error would affect both systems the same (and hence 
"correction" would not be useful) , but survey error in locating, 
say the scoring radar, as well as scoring radar tracking error, 
would inflate calculated CEP (and thus correction for these errors 
should be made, if detected, by changing p appropriately) . 

There is yet another source of error affecting calculated 
CEP's: that involved in estimating the MPI, y. These errors 

would be due to scoring system errors, such as scoring radar 
tracking errors. At present it appears that adjusting calculated 
CEP so that it does not include such scoring system errors (i.e. 
"discounting" calculated CEP so the delivery system is not incor- 
rectly "charged" with scoring system random effect errors) requires 
further theoretical development. We assume the scoring system has 
accuracy characteristics, when properly located near the drop point, 
such that serious error in calculated CEP does not result from ignor- 
ing these errors. 

In summary, using our radar bombscoring approach would 
involve determining inputs y and 2 to a CEP computation pro- 
cedure to be discussed in the next section. The vector y is 
estimated using output from the scoring system as discussed above. 

The matrix 2 is the sum of known variance-covariance components 
due to ballistic error, and delivery system random effects (if any) . 

4 . Computation of CEP Measures . 

Mathematically, our problem is to compute the median of 
the random variable L'L, given L has known bivariate normal 
distribution, L~N(y,2). This is a well known "damage" or "coverage" 
problem, and a large literature exists in connection with it. (See, 
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for example, the reference list of Johnson and Katz [6].) This 
problem amounts to finding a radius such that integration of the 
general bivariate normal distribution over an offset circle of 
that radius gives value 1/2. Computer routines for numerically 
approximating such integrals have been published, for example, by 
DiDonato and Jarnagin [3] . Tables of radius values for selected 
VI and ^ have also been published [4] . These tables, or an 
extension of them, might prove adequate for the CEP calculation 
problem. In what follows, we discuss an alternate approach. We 
consider first a special case (the non-central x^) then the more 
general case. 

If Y. were of the form (the offset circular normal 

situation), then as mentioned above, calculating CEP amounts to 
finding (or approximating) the median of a non-central distri- 

bution. Again, tables are available [5] or may be developed giving 
median values as a function of and X or, perhaps more 

naturally, as a function of cr^ and iJ • Alternatively, approxi- 
mate analytical expressions have been developed (see Johnson and 
Katz [6] , for example) ; a simple approximate median is given by 

CEP » 1.1774 0 (1+. 5000 X +.0817 X^) . (2) 

This expression can be viewed as a modification of the CEP formula 
(1) used for the non-offset circular normal case. Indeed, with 
X = 0 the offset circular normal situation specializes to the 
circular normal; formula (2) specializes to formula (1) as should 
be the case. 
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It is not difficult to develop a numerical procedure for 

this two degrees of freedom case which is simple enough to be 

implemented with a modestly sophisticated calculator. It is well 

known [5] that the non-central distribution can be expressed 

as a mixture of central distributions. Let denote the 

cumulative distribution function (CDF) for the x^ distribution 

with V degrees of freedom. Then the x.^ -.v CDF is given by 

( 2 ^ A ) 



F(z) 



00 



I 

k=0 



-X. k 

~lTl ^2+2k * 



(3) 



Repeated application of integration by parts to the integral 
representation of ®2+2k Yields the representation 






f)' 



(4) 



Replacing G in (3) by this representation, one obtains 



oo X 

F(z) =1-1 I 
k=0 j=0 



,k -X 
X e 

kJ 






(5) 



The computation problem can now be stated as follows; find n > 0 
such that 



-x-n 7 ? ni 1 

in 1^! jio 2 



k=0 



or (changing order of summation) , 

^"'X-n 



oo j j-1 k 



V ^ \ = I 

L ^ k ' 2 * 

j=l k=0 



CEP is then given by 



( 6 ) 



(7) 



CEP = a/2n 



( 8 ) 
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A ssairch iroutin© can ©asily b© giv©n to find n for any giv©n 
X > 0, and cons©qu©ntly to giv© CEP for any a > 0, using (8). 
Such a siinpl© expansion can b© obtained for any non~c©ntral chi~ 

square distribution with even degrees of freedom; otherwise the 
expansion would be such that calculation with a computer might 

be required. 



For the non-circular case, L'L cannot be simply normalized 
to a x'^ random variable. We now propose an alternative to the 
computer solutions mentioned at the beginning of this section. 

The approach is to approximate the distribution of L'L with that 
of a function of a certain non-central chi-square distributed 
random variable. This procedure is similar to that suggested by 
Pearson [10] for approximating non-central chi-square random vari- 
able in terms of functions of central random variables. In 

it, parameters in the approximating function are fitted by moments. 
Let us suppose, then, that 






(9) 



where means "approximately distributed as," the "scale factor" 

c and "degrees of freedom" v are parameters whose values must 
be determined, and where the non-centrality X is of the form 
X = y ' y . This form for non-centrality is suggested by the 
formula for the offset circular normal case. 



Assume L~N 



ff^ll 




'Oi 0 'I ' 




/ 


[0 a^JJ 



so L'L = + R 



and 



DVan" X'f2 






and R^/a^ ^ X/^ 






It follows that 



E(L'L) = a^E(D^/a^) + o^E(R^/a^) 



= a^(l+yj/2a^) + ' 



(10) 



T C 



and 



V(L'L) = aJv(DVo^) + 

= 2aJ[(l+y2/a2) + 2a^(l+y2/a2) , dD 

On the other hand, by assumption (9) it follows that 

E(L'L) i c(v+X) 

= c:(v + ^ (y^+y^)) (12) 

V(L'L) i 2C^ (V+2X) 

= 2C^ (V + i (y^+y^) ) • (13) 

Equating expression (12) with (10) , and (13) with (11) , and 
solving simultaneously for c and v, we obtain 

aj (aj+yj) + o^{o^+\i^) 

C = — — j j -- - (14) 

(o]^+y]^) + (‘^2'*'^2) 

(a^+y^) + (a^+y^) 

V = 1 + -^!^ — = — — - — — (15) 

a|(a|+y|) + 02 (O 2 +y 2 ) 

In practice, for associated CEP computation, it might be 
more tractible to round the value of v computed in (15) to the 
nearest integer value v*, and to take the corresponding c 
value, c = (0^+a|)/v*. For a number of example y and X values 
we examined, v* ranged between 1 and 3, which appears quite 
reasonable in comparison with the offset circular case. Indeed, 
acceptable accuracy might be achieved in many applications by 
taking v = 2 and c = (0^+0|)/2. This choice would allow use 
of the simple calculation approach of equations (6) and (7). Note 
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that in the offset circular case (i.e. a^=a|). 



use of formulas 



(14) and (15) would give 



earlier: With c = 



is exactly 



I 7 . 



(2,y *y/2c) * 



the correct distribution of L'L 
and V = 2, the distribution of 



given 

L'L/C 



5 . Summary . 

The approach we have described consists of using radar 
track, and other information, depending on the systems used, to 
estimate the MPI which would result if bombs were actually released. 
This vector, u, together with known ballistic and delivery system 
variance/covariance, S, can then be used in estimating the CEP 
associated with these values. The computation is based on approxi- 
mating the distribution of squared radial miss distance with a 
non-central chi-square distribution whose parameters are functions 
of the estimated y and known S. Several approaches to computa- 
tion of medians of such non-central chi-square distributions have 
been mentioned, and one procedure for use with the "even degrees 
of freedom" case was presented. 

The method of estimation of the MPI for a given bomb run 
depends upon the systems being used for delivery and scoring. 

The combinations we discussed are indicated in the following dia- 
gram with the cells marked "x." 

System Configuration Computation Method 



Delivery 


Scoring 


Ballistic Computation 


Bombing Tables 


Pilot controlled 


Separate 


X 


X 


System controlled 


Separate 


X 




System controlled 


Common 


X 
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The roles of various potential error sources may be different 
in the 'Separate " scoring system case from that in the "common" 
system. Example fixed, random and unaccounted effects are shown 
below. 





Fixed Effect 


Random Effect 


Unaccounted 


Common deliv- 
ery and 
scoring 
system 


Pilot error 

Delivery system 
error (such as 
not dropping 
at proper time 

System inability 
to achieve 
exactly desire 
drop conditior 


Ballistic error 

0 

7 

id 

i 


Air density profile 
Drag curve 
Radar survey error 
Radar track error 


Separate scor- 
ing system 


Survey error 

Delivery 

System 

Tracking 

Error 


Scoring system 
tracking error 


Errors in drop 
condition inputs 
(winds, drag, etc.) 



Several problems remain with respect to details of 
implementing the approach suggested herein, and effort should 
be directed toward their solution. First, what is the best way 
to calculate (or approximate) the median of L'L in the non- 
circular offset normal case, and what accuracies can be expected 
with each candidate method? Second, how should calculated CEP's 
be "discounted" to remove errors and variance components which 
are due to the scoring system, rather than the delivery system? 
Third, what is the best way to combine CEP estimates made in 
individual runs into a single CEP measure for a "sortie?" We 
plan to consider these and related questions in future efforts. 
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